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Towards  Shorter  Wavelength  Lasers  and 
Breaking  the  1000  A Barrier- I II 
( Recombination) 


1.  Introduction 

Discussions  on  the  possibility  of  VUV  and  X-Ray  lasers  can  be  found 
in  the  literature  as  long  as  ten  years  ago1*2.  The  interest  in  the  sub- 
ject, however,  was  enhanced  in  19 6?  when  the  possibility  of  photon  pump- 
ing as  a pos  ible  mechanism  for  the  realization  of  VUV  (X  < 1000)  and 
X-Ray  lasers3  (X  < 10  A)  was  discussed.  Meanwhile,  there  have  arisen 
a variety  of  proposals  and  schemes  on  this  subject.  The  electron  and 
photon  pumping  mechanisms  were  discussed  in  Ref.  (4)  and  the  Pound- 
continuum  emission  as  a source  for  short  wave  lasers  was  discussed  in 
Ref.  (5).  In  this  report  the  recombining  plasma  is  considered  as  a 
medium  for  the  amplification  of  short  wave  lasers  (X  < 1000  A). 

Gudzenko  and  Shelepin0  were  first  to  point  out  the  possibility  of 
radiation  amplification  in  a recombining  plasma.  Their  calculations  for 
an  optically  thin  hydrogen  plasma  based  on  non- equilibrium  plasma  con- 
ditions showed  inversion  and  reasonable  gain  in  the  visible.  Their 
method  basically  follows  the  standard  work  of  Bates,  et  al.7,  who  have 
calculated  the  steady  state  values  of  the  excited  state  densities  for 
optically  thin  and  thick  hydrogen  and  hydrogen  like  plasmas.  These 

*7 

calculations  were  carried  out  for  given  electron  densities  and  tempera- 
tures with  the  motivation  to  obtain  the  total  recombination  coefficient 
and  the  deviation  of  the  excited  states  from  their  equilibrium  (Saha) 
densities.  An  optically  thick  calculation  for  laser  purposes  was  also 
carried  out  by  Gordiets  et  al.s,  for  a hydrogenic  plasma.  Bohn5, 
recently,  have  calculated  inversion  densities  and  gain  coefficients  for 
hydrogen  like  ions  with  emphasis  on  short  wave  lasers.  His  calculations9 
Note:  Manuscript  submitted  February  14,  1975. 
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follow  the  standard  Bates  et  al.7  method.  A similar  calculation  for 
hydrogen  plasma  is  carried  out  by  Skorupski  and  Suckewer*0  with  emphasis 
on  the  effect  of  heavy  particle  (hydrogen  atom)  collisions  on  the  ex- 
cited state  populations. 

For  hydrogen-like  ions  the  energy  spacing  between  discrete  levels 
increases  as  Z2  (where  Z is  the  nuclear  charge).  Therefore,  shorter 
wavelength  lasers  (X  < 1000  A)  can  be  realized  within  a few  times  ionized 
hydrogen- like  plasmas.  For  example  the  n = 3 to  n = 2 transition  in 
hydrogen  (where  n and  m are  the  quantum  numbers)  corresponds  to  6^62  A, 
which  is  the  well  known  H line.  In  He  II  the  same  transition  is  1640  A 
while  in  Li  III  it  corresponds  to  729  A. 

A time  dependent  solution  for  the  inversion  density  and  the  gain 
coefficient  with  a self  consistent  treatment  of  radiation,  we  believe, 
is  ? step  forward  in  the  recombination  plasma  model  calculations.  This 
has  been  carried  out  by  Jones  and  Ali11.  It  should  be  pointed  out, 
however,  that  a time  dependent  calculation  for  He  IT.  with  emphasis  on 
1640  A line  has  been  carried  out  by  Gudzenko,  et  al.12  for  a helium 
plasma  where  the  effect  of  charge  exchange  between  an  excited  Helium  ion 
state  and  a Helium  atom  is  considered.  They  have  also  considered13  some 
aspects  of  cw  operation  of  a recombining  plasma.  Inversion  between  the 
resonance  and  the  ground-state  of  hydrogen  in  a two  level  atom  approxima- 
tion has  been  considered  analytically  by  Peyrand  and  Peyrand14. 

One  important  aspect  of  the  recombination  plasma  as  a source  of 
radiation  amplification  is  that  the  electron  cooling  must  be  faster  than 
the  plasma  recombination.  Some  aspects  of  this  has  been  considered  by 
Pert  and  Ramsden15.  On  the  other  hand,  Irons  and  Peacock16  indicate 

+5 

the  presence  of  some  population  inversion  in  C ion  in  an  expanding 
laser  produced  plasma. 

In  this  report,  the  physics  of  recombination  as  a source  for  short 
wave  lasers  (X  < 1000  A)  is  considered  and  some  results  of  a model 
calculations  are  presented. 
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2.0  The  Basic  Physics  of  Recombination  Plasmas 


The  process  of  recombination  which  signifies  the  relaxation  phenom- 
enon in  an  ionized  medium  has  been  of  great  interest  to  astrophysicist 
and  plasma  spectroscopists.  Radiation,  as  line  spectra  and  continuum, 
accompany  the  recombination  process.  A quantitative  measurement  of  these 
spectra  enables  one  to  describe  the  physical  state  of  a plasma. 

The  astrophysicists  have  generally  investigated  plasmas  with  a low 
degree  of  ionization,  where  collisional  processes  (N£  « 10s  cm-3)  are 
not  significant.  Radiative  capture  and  cascade  are  the  dominant  pro- 
cesses which  could  describe  the  line  intensities  from  gaseous  nebulae, 
for  example.  Baker  and  Menzel17*18  originally  approached  this  problem 
by  calculating  the  line  intensities  from  a gaseous  nebulae  for  two  well 
known  cases:  plasmas  1)  optically  thin  and  2)  optically  thick  to  Lyman 
series  radiation.  Application  of  their  techniques  to  astrophysical 
problems  with  further  refinements  have  been  carried  out  by  numerous  in- 
vestigators: 1)  Burgess19,  taking  explicit  account  of  the  orbital  angular 
momentum  and  2)  Seaton20  taking  into  consideration  the  three-body  recom- 
bination for  higher  quantum  levels. 

It  is  of  interest  at  this  juncture  to  point  out  that  based  on  purely 

radiative  recombination  and  cascade  (N  « 108  cm-3)  one  can  not  obtain 

e 

inversion  in  a hydrogen  or  hydrogen-like  plasma.  This  is  because  direct 

recombination  favors  lower  quantum  levels.  For  example,  for  hydrogen  the 

recombination  coefficient  CC  (n)  to  level  n is21 

r ' 

a (n)  cel  l , (1) 

n NT 

when  the  mean  thermal  energy  of  electrons,  kT,  is  small  compared  to  the 
ionization  energy.  However,  when  the  thermal  energy  is  larger  than  the 
ionization  energy,  then21 


a (n) 
r 


(2) 
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again  favoring  the  lower  levels.  Furthermore,  detailed  calculations  by 
Drawin  and  Emard  for  the  instantaneous  population  densities  of  ex- 
cited states  show  lack  of  inversion.  For  example,  for  T — 4000°K  and 
Ne  - 1013  cm  3 the  departure  coefficient  from  Saha-Equilibrium  for 
n = 4 and  n = 3 are  6.5  X 10  2 and  8.25  X 10  2 respectively,  indicating 
that  no  inversion  is  possible.  It  should  be  pointed  out,  however,  that 
even  if  inversion  occurs  at  such  a low  density,  gain  is  insufficient  for 
laser  action. 

For  dense  plasmas,  however,  with  higher  degrees  of  ionization,  the 
recombination  is  essentially  collisional  and  three-body  recombination 
dominates.  The  pioneering  work  of  Rates  et  al.7  with  their  extension 
of  recombination  to  plasmas  with  higher  electron  densities  and  thereby 
including  three-body  processes,  have  been  followed  by  numerous  investi- 
gators with  interest  in  plasma  physics  and  diagnostics  of  laboratory 
plasmas.  These  calculations  have  either  solved  a large  set  of  rate 
equations7’23-25  which  describe  the  population  densities  of  excited 
states,  or  are  based  on  a variational  model26’27  which  yields  the  total 
recombination  coefficients.  Therefore,  for  laser  purposes,  with  a 
recombining,  high  density  plasma,  the  method  of  Bates  et  al.7  is  essen- 
tially the  starting  point,  since  it  can  solve  for  the  population 
densities  of  the  excited  states  and  their  method  has  been  followed  by 
many6’9  »10 . 

On  the  other  hand,  the  three-body  recombination  coefficient,  CC  (n), 

3 

has  a dependence  on  the  quantum  number  favoring  the  higher  levels.  This 

coefficient  can  be  obtained  from  the  ionization  rate  coefficient,  S , 

n 

through  the  detailed  balance  (see  reaction  3) 


N + N + NZ  **  NZ_1{n)  + N , 
e e ' ' e 


where  Ng  is  the  electron  density,  N is  the  ion  density,  Z times  ionized, 
and  NZ_1(n)  is  the  population  density  of  level  n within  the  ion  which 


is  Z - 1 times  ionized. 


The  ionization  rate  coefficient  can  be  expressed  as 


S 

n 


8.0  X 10" V 

exp 

ZaEu  \!  T 
H 1 o 


where  is  the  ionization  energy  of  hydrogen.  Here  Tg  and  are  in 
units  of  eV.  This  expression  is  obtained  using  the  excitation  rate 
coefficient4  given  by  Seaton  with  an  oscillator  strength  of  unity  and  an 
average  gaunt  factor  of  j . Bates,  et  al.  using  an  expression  of  this 
kind  indicate  that  (4)  is  reasonable  within  a factor  of  2 for  helium. 
Using  detailed  balance  one  obtains 

N NZ  S 

e _ n 

NZ-1(n)  a (n) 

3 

However,  the  left  hand  side  of  (5)  is  the  well  known  Saha  distribution28 
which  can  be  expressed  as 


(5) 


N 

e 

N^U) 


3.02  X 1021  T 


3/a  J *%\ 


exp  l - 


n2T 


(6) 


where  gZ  is  the  statistical  weight  of  the  ground  state  of  NZ.  Using  (4)  - 
(6)  one  obtains 


a (n)  = 2.65  X 10' 

3 


-27 


(cm6/ 


sec : 


(7) 


It  is  obvious  that  the  three  body  recombination  rate  coefficient  for  a 
given  Z is  proportional  to  n4  and  the  lower  the  temperature,  the  higher 
the  recombination  rate.  Furthermore,  for  the  same  quantum  level,  the 
recombination  rate  coefficient  is  inversely  proportional  to  Z . 
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2.1  Hydrogen-like  Energy  Level  and  Gain 


The  wavelength  for  a transition  between  two  hydrogen-like  energy 
levels  n and  m is 


X 


nm 


1.24  X 10" 4 n2m2 
Z2E^  n2  - ir2 


(8) 


where  E„  is  in  units  of  eV.  Consider,  for  example,  a 3 2 

H o o 

transition.  This  transition  is  6562  A in  hydrogen,  l64o  A in  He  II  and 
729  A in  Li  III.  Thus  the  breaking  of  the  1000  A barrier  can  be  achieved 
in  a recombining  plasma  where  Li  is  three  times  ionized.  Shorter  wave- 
length lasers  can  be  obtained  if  one  considers  transitions  with  An  > 1 
in  the  above  ions. 

The  inversion  density  between  two  discrete  levels  u and  Z with 

Nu 

densities  N and  N„,  respectively  is , where  g and  g.  are  the 

U Xf  U Yj 

8U  -1  4 

corresponding  statistical  weights.  The  gain  coefficient,  in  cm  is 


a = 10  12  A „ 

u Z 


(9) 


where  is  the  spontaneous  transition  rate  and  AXu^  is  the  width  of  the 
transition.  One  can  express  the  gain  in  Eq.  (9)  for  a hydrogen- like  ion 
in  terms  of  the  quantum  numbers  of  the  levels,  the  oscillator  strength 
of  the  transition  and  the  nuclear  charge.  This  can  be  done  by  utilizing 
the  expression28  for  the  transition  rate  i.e. 


A 


u Z 


•7  ® Z 

= 4.3  x 107  — f ,z 


uf  “II 


JL_ 

r 


(10) 


Thus,  using  (8)  and  ( 10)  into  (9)  one  obtains 
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N - N.  — ^ 

“ i »i  a z* 


u2  (u2  - J J2)2 


In  most  cases  of  interest,  AX  «,  is  Stark  broadened  and  must  be  calculated 

u£ 

for  each  electron  density.  However,  when  AXu^  is  Doppler  broadened  it  is 
proportional  to  \ . and  thus  to  — . Therefore,  the  gain  in  terms  of 

£-2  J 

quantum  levels  as  given  in  (ll)  would  be  proportional  to — . Using  the 

Z2  _ 

generally  accepted  scaling  for  hydrogen-like  ions  i.e.  N = n Z and 

G 6 

= Z 9,  where  ng  and  8g  are  the  electron  density  and  temperature  for 
hydrogen,  one  obtains  the  scaled  gain  to  be  proportional  to  Z7.  This 
scaling  for  a Doppler  broadened  line  and  T^  = Tg  is 


' °B  ' 2 


, M 

7 2 


% 'H, 


2.2  Inversion  and  Recombination  Relaxation  Time 

When  a plasma  is  completely  ionized  i.e.,  when  the  constituents  are 
completely  stripped  ions  and  electrons,  the  plasma  will  relax  unless  the 
electrons  are  continuously  heated  by  an  external  source.  After  a short 
time  the  plasma  reaches  a quasi  equilibrium  state.  During  this  early 
stage  of  recombination  one  is  expected  to  see  laser  action  due  to  the 
filling  of  the  higher  states.  This  is  because,  the  three  body  recombina- 
tion clearly  favors  the  higher  quantum  levels  (see  Eq.  7)»  To  calculate 
the  total  recombination  rate  it  is  instructive  to  point  out  that  there  is 
a quantum  level23 , n*,  above  which  all  discrete  excited  levels  are  in 
equilibrium  with  the  free  electrons.  Furthermore29,  the  number  of 
transitions  transferred  from  continuum  downward  thru  this  level,  equals 
the  number  of  transitions  upward  from  below  the  same  level.  This  quantum 
state  is  given  by 
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and  can  be  arrived  at  by  the  following  argument.  Since  transitions  among 
the  nearest  neighbors  are  the  most  important  transitions,  one  obtains 

N , Y , , -NX 
n+i  n+i,n  n-i  n-i,n 


where  X /is  the  excitation  rate  from  n to  n;  and  Y / is  the  corre- 
n,n  n,n 

sponding  de-excitation  rate.  Equation  (13)  implies  that 


or  E - T 
n e 

sum  of  Eq . 

(t^-1  » 2.65  X lCf28  (Z2EHf /2(Te)-4-V  , (15) 


and  n* 


N 


r2v- 


"h 


Thus.  the  total  recombination  rate  is  the 


(7)  from  n = 1 to  n = n4:"  which  results  in 


where  = 2 was  utilized.  This  result  differs  (for  Z = 1)  from  corre- 
sponding results  in  Refs.  29  and  50  by  a factor  of  2.35  and  1.5,  re- 
spectively. These  differences  arise  primarily  because  of  the  choice  of 
the  cross  sections  used.  However,  the  more  rigorous  result  derived  by 
Pitaevskii  and  by  Gurevich  and  Pitaevskii  (see  Ref.  30)  is 


aT  = ^rr(2n)1/'2 


e10Z3 

'^(kT)9/2 


log  A 


which  reduces,  with  log  A * 1,  to 


(16a) 


0T  = 8.75  X IQ"57  Z3 

3 (T  )4*5 

e' 


( cm6/ sec) 


(16b) 


The  range  of  applicability  of  this  relation  is  for  Z2E„  » T or  n*  » 1. 

n 6 

For  hydrogen  the  validity  of  this  relation  is  for  T^  ^ 0.26  eV. 
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However,  a total  three  body  recombination  rate  coefficient  will  have 
an  entirely  different  temperature  dependence  and  indeed  even  a different 
value  if  one  utilizes  the  ionization  rate  coefficients  reported  by 
McWhirter31 . Consider  the  following  expression  for  the  ionization 
coefficient  of  species  Z 

(Te/EZ_1)1/2 


1,  with  ionization  energy  E2"1  = Z2E^. 


,z-i 


9 X 10" 


,2-r 


(ez-1)3/s 


exp  ( 


(17a) 


4.88  + 


V TC 


,Z-1 


or 


.z-l 


1,1  X 10" 


/ 2 


(E2"1) 


*3/2 


exp 


(17b) 


o + ■ 


,z-l 


These  expressions  are  valid"31  over  a wide  range  i.e.,  .02  < 


,z-l 


< 100, 


which  clearly  encompasses  the  region  of  the  validity  of(l6b).  The 
corresponding  three  body  recombination  rate  coefficient,  using  (17a)  is 


T 

a = 

3 


1.5  x 10"27  (t  /e2"1)172 


(E  T ) 
e 


3 /p 


(18) 


4.88  + 


,Z~1 


which  for  T /E 
e 


z-l 


« 1 clearly  has  an  entirely  different  temperature  and 


Z dependence,  that  is  C!^  « Z ° Tg  . For  a comparison  between  ( 18)  and 
(l6b) let  T = 1 eV,  Z = 2.  The  rate  coefficients  will  be  10  27  cm6/sec 
and  70  X 10~27 


cm^/sec  respectively.  Tliis  discrepancy, which  is  even 
larger  than  two  orders  of  magnitude  in  certain  instances  will  have  a 
profound  result  on  resonance  line  transitions  if  one  considers  how  fast 
the  hydrogenic  ion  becomes  a helium- like  ion. 

Starting  from  completely  stripped  ions,  the  early  stage  of  the 
recombination  results  in  inversion  and  hence  laser  action  among  the 
highly  excited  and  discrete  states  of  the  hydrogen-like  ions  (see  Fig.  1 
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for  inversion  densities  in  hydrogen).  Laser  action  generally  terminates 
after  times  short  compared  to  the  life  times  of  the  transitions.  These 
times  are  of  the  order  of  the  electron  de-excitation  times  of  the  upper 
laser  levels  or  times  when  equilibrium  conditions  are  attained.  The  time 
required  for  level  n(n  r 1)  co  be  in  equilibrium  with  the  higher  levels 
n7  > n and  the  free  electrons  is  equal  to  the  inverse  of  the  sum  of  the 
excitation  rates  to  all  higher  levels,  or  approximately,  to  the  next 
higher  level.  For  a hydrogen- like  ion,  the  excitation  rate  coefficient 
from  level  n to  n7  is’1 


X / 
nn 


l.o  X 10"~  f t 
nn 

AEnn  ' \|  ? 
nn  ] e 


(19a) 


where  T and  t E / are  in  eV.  Assuming  a unit  Gaunt  factor  reduces  and 
e nn 


2 Z^E. 


nn 


using  , and  T f / ^ 0.5  n (see  Ref.  26),  (19a)  becomes 

h x 10"' ‘ n4  / 2 Z2E„ 


X / = — 
nn  ,2 


exp 


Z E..  T 
H e 


n3T 


(19b) 


Thus,tne  time  for  level  n to  reach  equilibrium  with  the  upper  levels  is 


2.5  X 10= 


n4N 


Z%'|  Te  exp  ( 


ifV, 

* 3frt 

\ n T 


(20) 


This  expression  (20),  should  be  used  when,  in  general,  there  are  many 
discrete  levels  above  the  level  of  interest.  On  the  other  hand  when  only 
a few  levels  are  discrete  (see  the  discussion  on  the  reduction  of  tha 
ionization  limits  in  section  2.5)  then  the  de-excitation  time  from  one 
level  to  another  will  determine  the  termination  of  the  laser  action. 
Furthermore,  for  equilibration  down  to  the  ground  state  of  the  ion,  the 
de-excitation  rate  of  the  resonance  level  determines  the  required  time 
for  LTE  to  set  in.  In  this  case  the  time  of  interest  is 
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T 

21 


1 


(21) 


where  Y , is  the  electron  impact  de-excitation  rate  coefficient  from 
n n 

level  n'  to  n(n/  > n)  and  it  can  be  obtained  from  (19)  via  detailed 
balance. 

2.3  Excited  Level  Population 

The  laser  action  in  a recombining  plasma  occurs  between  a pair  of 
excited  states  of  an  atom  or  an  atomic  ion.  It  can  also  take  place 
between  an  excited  state  and  .ae  ground  state  of  the  atom  or  the  atomic 
ion.  Therefore,  a knowledge  of  the  population  densities  of  these  states 
is  essential  for  the  laser  power  density  calculations,  the  requirements 
lor  an  efficient  laser  and  the  optimum  conditions  for  lasing.  The 
population  density  of  a level  can  be  obtained,  in  general,  either  by 
following  the  techniques  of  Bates  et  al.7  which  gives  the  steady  state 
solutions  for  a set  of  rate  equations, or  by  solving  the  rate  equations 
as  a function  of  time.  The  steady  state  solution  ignores  the  early  time 
spikes  of  lasers  which  may  occur  as  the  levels  are  being  filled  before 
they  eventually  attain  a steady  state.  Therefore,  the  steady  state 
underestimates  the  rate  at  which  the  charge  disappears  at  early  times 
where  the  levels  are  being  filled.  Furthermore,  a steady  state  solution 
for  a generalized  hydrogenic  system  can  not  account  for  the  loss  of  the 
N stage  of  ionization.  However,  a time  dependent  model  which  also 
solves11  for  the  radiation  densities  connecting  any  two  pairs  of  levels 
will  self  consistently  describe  the  actual  time  development  of  the  levels. 
In  fact  a model  which  also  considers  the  time  variation  of  the  electron 
temperature  due  to  the  inelastic  collisions  and  the  dynamic  expansion  of 
the  plasma  is  much  superior  to  the  steady  state  solution.  The  processes 
affecting  the  population  density  of  an  exciced  level  are  three-body 
recombination,  radiative  recombination,  spontaneous  decay  and  stimulated 
emission  from  higher  levels,  electron  superelastic  de-excitation  of 
higher  levels,  radiative  absorption  from  lower  levels  and  electron  impact 
excitation  from  lower  levels.  All  the  above  processes  increase  the 


11 


population  density  of  a given  level.  In  contrast,  the  processes  which 
decrease  the  population  density  of  a level  are:  radiative  decay,  stimu- 
lated emission,  electron  impact  de-excitation  to  lower  levels,  electron 
impact  excitation  which  raise  the  electron  to  higher  levels,  ionization 
by  electron  impact,  photoabsorption  and  ionization. 

Thus, in  principle  one  must  solve  for  a large  number  of  excited 
states.  However,  the  number  of  states  is  limited  because  of  the  re- 
duction in  the  ionization  limit  in  a plasma.  The  last  discrete  level, 

. r>  c. 


v 1S‘  - 
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\ AE 


z-i 


1/2 


(22) 


where  AE 


z-i 


is  the  reduction  in  the  ionization  limit  and  is  given  by 


AE 


z-i 


(25) 
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In  (A5'  tp  is  the  Debye  radius  which  is  defined  as 


°D  = 6,9  ' ~ 

\ N / 


1 /£ 


(2M 


where  T is  the  electron  temperature  in  degrees  Kelvin  and^* 


N7  = N + I Z2NZ  . (25) 

e z 

In  general,  however,  most  plasmas  are  not  completely  ionized.  There- 
fore, the  population  densities  of  the  excited  states  of  interest  may  be 
affected  by  other  processes  which  must  be  included.  Among  these;  is 
collisions  of  heavy  particles  (atoms  or  ions)  with  the  excited  state. 

These  collisions  may  excite,  de-excite  and  cause  charge  transfer  to  or 
from  the  level  of  interest. 
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2 .h  The  Electron  Temperature  and  the  Average  Charge  State  in  a Plasma 


A knowledge  of  the  electron  temperature  and  its  evolution  is  of 
considerable  interest  in  the  recombination  calculations.  For  example  the 
cooling  of  the  electrons  increase  the  three-body  recombination  rate  to  a 


given  level,  at  a faster  rate 


— — 'N  than  the  increase  in  the  level  de- 
T2  / 

a ' 


excitation  rate  by  electrons  y-  Furthermore,  for  high  Z ions,  the 

cooling  depletes  the  ground  stateGat  a faster  rate  f— - — \ ; thereby  in- 

\T4*5/ 

creasing  the  inversion  density  and  the  duration  for  n = 2 to  1 transitions. 


Therefore  the  time  history  of  the  electron  temperature  should  be  known 
for  the  plasma  of  interest.  In  this  section  a semi-quantitativa  approach 
is  taken  to  estimate  the  electron  temperature  and  the  charge  state  in  a 
plasma  in  relation  to  the  energy  input,  especially  for  the  laser  produced 
plasmas.  High  temperature  and  high  density  plasmas  are  obtained  in 

33  34 

several  ways.  The  most  notable  and  novel  is  the  laser  produced  plasma  9 
For  a review  article  on  this  subject  see  Ref.  35.  One  of  the  important 
aspects  of  this  novel  approach  is  that  high  electron  densities  n£  » 

1010  cm  3 and  high  electron  temperatures  T =“  keV  can  be  achieved.  How- 
ever, plasmas  with  high  electron  densities  can  also  be  obtained  through 
shock  tubes30  3 &,  coaxial  accelerators39  and  by  high  energy  electron 
beams40.  It  is  interesting  to  point  out  at  this  juncture  that  the  state 
of  the  plasma  in  shock  tubes  e.g., whether  LTE  is  established  or  not  have 
been  a center  of  controversy  among  the  investigators3 ° 39 . It  would  be 
desirable  to  see  if  population  inversion  occurs  in  these  transient 
plasmas  and  whether  or  not  this  mechanism  may  be  the  reason  for  the  dis- 
agreement^0 39 . The  authors  plan  to  review  this  in  the  near  future. 

In  a laser  produced  plasma,  the  electrons  gain  energy  from  the 
incident  laser  beam  through  inverse  bremstrahlung  (in  a collisional 
plasma)  and  cool  through  ionization,  excitation,  electron-electron  and 
electron- ion  collisions.  The  three  body  recombination  heats  the  electrons 
in  a recombining  plasma,  however,  as  the  plasma  expands  into  a vacuum, 
the  electrons  cool  adiabatically . Thus,  the  electron  temperature  can  be 
described  provided  all  relevant  processes  are  considered.  The  energy 


15 


density,  w,  deposited  in  a target  with  charge  Z to  produce  Ng  electrons 
with  temperature  T^  can  be  written  as 

/ » 

w = - N T N +V'  N E 

2 e a 2 1 / z-x  z-x  z-x 

Z=  1 2=1 

(26) 

/ 

+ > N I , 

/ i z-x  z-x 

z=x 


where  T.  is  the  average  ion  temperature  (in  eV),  E is  the  excitation 
t Z“  x 

energy  of  the  internal  states  of  species,  Z - 1 and  I is  the  ioniza- 
tion energy  of  the  same  species.  The  prime  on  the  summations  indicates 
the  last  stage  of  ionization  desired.  For  example,  in  a completely  stripped 
helium,  the  summation  is  up  to  Z = 2.  Thus,  for  the  case  of  completely 
stripped  ions,  the  energy  density  is 

2 

E + 

z-x  ^ ZT 

e z=l  .3 

This  coupled  with  the  degree  of  ionization  and  average  charge  state 
calculations  can,  in  principle,  give  the  electron  temperature  phenom- 
enologically. On  the  other  hand,  the  temperature  evolution  can  be 
calculated  from  the  first  electron  stripping  to  the  last  stages  of 
ionization  in  a laser  produced  plasma  from  the  laser  absorption  principle. 

The  average  charge  state  in  a plasma  can  be  calculated  following 
Raizer30 »‘u  by  assuming  that  the  ion  densities  and  the  ionization 
potentials  are  continuous  functions  of  the  ionic  charge  Z.  This  can 
be  expressed  as 


w = — N T 
^ e e 


ZT 


3 ZT 


V 

C-J 

z=  l 


N = N + — 
2 "-1  dZ 


) 


or 
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= 1 + 


1 dN 


N , dZ 
z-1 


However,  — — is  zero  at  the  peak  of  the  ionization.  Therefore,  for  a 
dZ 

collision  dominated  plasma  expression  (28)  can  be  related  to  Saha 
equation  (see  Eq.  6 for  n = l)  i.e. 


N = 6 X Itf21  T 


where  T is  in  units  of  eV,  the  statistical  weights  are  set  to  unity, 
and  I is  the  average  ionization  energy.  Thus  for  a given  electron  density 
and  temperature  one  finds  the  average  ionization  energy  and  thus  the 
average  charge  state  using  values  of  the  ionization  potentials  . The 
charge  state,  Z,  thus  obtained  should  correspond  to  the  average  charge 

Z by  simply  adding30’41  one  half  to  Z. 

As  an  example,  consider  carbon  ionized  six  times  which  implies  that 

f s = lj-90  eV  and  assume  an  electron  density  of  10s2  cm  . Using  (29) 

one  obtains  T « 80  eV.  This  can  be  coupled  with  (27)  to  give  the 
e . 

required  or  the  minimum  energy  input.  Assuming  that  Tg  T\,  the  energy 

density  can  be  approximated  for  the  above  condition  by 

6 

w ^ — N T 1.2  + V I 

2 e e Lm4  z~± 


1.2  + 


e z=  l 


w - 6 X 105  J/cm3  . 


This  energy  should  then  be  delivered  in  times  short  compared  to  the 
recombination  time  and  for  lasers  in  the  resonance  state,  the  energy 
should  be  delivered  in  times  short  compared  to  the  electron  de-excitation 
time  of  the  resonance  level.  For  carbon,  this  time  is  ~ 10  sec  l see 
Eq.  21).  This  implies  that  6 x 1017  watts/cm3  is  needed  for  one  cubic 
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centimeter.  However,  if  one  considers  a carbon  sample  of  1 cm  long  and 
1 u wide  then  the  power  required  will  be  6 x 1C?  watts. 

The  assumption  that  the  electron  temperature  is  equal  to  the  ion 
temperature  over  estimates  the  energy  input.  The  equipartition  time 
between  electrons  and  ions  is?*~ 


I 

ei 


A 

5 x lo7  — 
z2n 


(50) 


where  A is  the  atomic  number  and  T is  in  eV.  For  T =•  80  eV,  Z = 6 
z e e 

and  !>  = lO^*"  cm  one  obtains  t . = ^ X 10  sec.  This  time  is  longer 

e gl 

than  the  duration  of  inversion. 


Quantitative  Estimates  of  Inversion  Density  and  Gain 

In  this  section  quantitative  estimates  of  inversion  density  and  gain 

for  resonance  line  transitions  in  hydrogen  like  ions  will  be  given. 

Consider  for  simplicity  a completely  stripped  carbon  with  Tg  = 100  eV 

and  N = lO2*  cm  3.  Using  equations  22-25  one  obtains  n,  =*  ^ which 
e d 

implies  that  only  few  discrete  levels  are  observable.  Furthermore,  the 
upper  level  of  the  resonance  line  i.e.  n - 2 will  be  in  equilibrium  with 
the  higher  discrete  levels  and  the  free  electrons.  This  can  be  obtained 
using  the  criterion  for  equilibration28  which  states  that 


7 X 1018  Zn 


N 5 
e 


/ T 
( 

V 


(51) 


However,  level  2 is  not  in  equilibrium  with  the  ground  state  and  in- 
version is  established  between  n = 2 and  n = 1 at  the  early  stage  of 
recombination  simple  because  the  recombination  rate  to  n = 2 is  larger 
than  the  corresponding  rate  to  n = 1.  In  addition,  level  n = 1 will  be 
depleted  towards  the  formation  of  the  helium-like  ion  and  the  establishment 
of  the  collisional  ionization  equilibrium. 

The  inversion  density  can  be  obtained  by  using  the  equilibrium  value 

for  n = 2.  Using  Eq.  (6)  one  obtains  N =“  4 x 1019  cm"3,  which  is 

2 

approximately  equal  to  the  inversion  density.  The  gain,  on  the  other  hand, 


1 6 


is  a 2:  28  cm”1,  where  equation  (9)  is  utilized  with  the  Stark  width  of 
the  transition.  The  full  Stark  width43  is 

n?)K2/*  • iJ  . (32) 

TT  H C 

Using  Eq.  (32)  one  obtains  AX  ^ 1 A,  which  is  much  larger  than  the 
Doppler  width  (assuming  that  T.  = Tg).  The  above  analytic  estimates 
agree  very  well  with  the  detailed11  code  calculations. 

However,  it  should  be  recognized  that  the  duration  of  inversion 
is  short  and  is  of  the  order  of  the  electron  de-excitation  rate  of 
level  2.  This  process  is  responsible  for  the  equilibration  of  level  1 
with  the  higher  levels  and  thus  will  terminate  any  laser  action. 

Figure  (2)  shows  the  inversion  density  and  the  gain  for  n - 2 to 
n = 1 transition  in  hydrogen- like  carbon  (Xs*  34  A)  for  several  initial 
electron  temperatures.  These  results  are  from  a detailed  model  calcula- 
tions 11 . However,  for  T^  100  eV,  the  simple  quantitative  analysis, 
given  in  this  section,  predicts  results  in  good  agreement  with  the  de- 
tailed calculation. 
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